Astronomy & Astrophysics manuscript no. 7962 


© ESO 2008 


February 2, 2008 





o 
o 

(N 
•*-> 
O 

O 
o 



OH" 

I 

o 



> 

in 



o 



X 



A synchrotron self-Compton model with low energy electron 
cut-off for the blazar S5 0716+714 

Olivia Tsang and J. G. Kirk 

Max-Planck-Institut-fiir Kernphysik, Saupfercheckweg 1, D-69117 Heidelberg, Germany 
Received . . . / Accepted . . . 

ABSTRACT 

Context. In a self-absorbed synchrotron source with power-law electrons, rapid inverse Compton cooling sets in when the 
brightness temperature of the source reaches Tb ~ 10 12 K. However, brightness temperatures inferred from observations of 
intra-day variable sources (IDV) are well above the "Compton catastrophe" limit. This can be understood if the underlying 
electron distribution cuts off at low energy. 

Aims. We examine the compatibility of the synchrotron and inverse Compton emission of an electron distribution with low- 
energy cut-off with that of IDV sources, using the observed spectral energy distribution of S5 0716+714 as an example. 
Methods. We compute the synchrotron self-Compton (SSC) spectrum of monoenergetic electrons and compare it to the observed 
spectral energy distribution (SED) of S5 0716+714. The hard radio spectrum is well-fitted by this model, and the optical data 
can be accommodated by a power-law extension to the electron spectrum. We therefore examine the scenario of an injection of 
electrons, which is a double power law in energy, with a hard low-energy component that does not contribute to the synchrotron 
opacity. 

Results. We show that the double power-law injection model is in good agreement with the observed SED of S5 0716+714. 
For intrinsic variability, we find that a Doppler factor of T> > 30 can explain the observed SED provided that low-frequency 
(< 32 GHz) emission originates from a larger region than the higher-frequency emission. To fit the entire spectrum, 23 > 65 is 
needed. We find the constraint imposed by induced Compton scattering at high Tb is insignificant in our model. 
Conclusions. We confirm that electron distribution with a low-energy cut-off can explain the high brightness temperature in 
compact radio sources. We show that synchrotron spectrum from such distributions naturally accounts for the observed hard 
radio continuum with a softer optical component, without the need for an inhomogeneous source. The required low energy 
electron distribution is compatible with a relativistic Maxwellian. 

Key words, galaxies: active - galaxies: high redshift - galaxies: jets - BL Lacertae objects: individual: S5 0716+714 



1. Introduction 

Observations of many extra-galactic radio sources have 
found rapid flux variations at radio frequency (e.g. 



Kedziora-Chudczer et al. 120011) . some of which fluctu- 
ate over a time scale of a day or less. They are re- 
ferred to as intra-day variable sources (IDV). The vari- 
ability time scale is often used to constrain the size of 
the source based on causality arguments. Using this con- 
straint, one can derive a v ariability brightness temperature 
(jWagner fc Witzellll995l) 



T mT = 4.5 x 10 1U F„ 



K 



(1) 



t ohs (l + z) 

where the flux density F Vl wavelength A, luminosity dis- 
tant g?l, and observed variability time scale t h s are mea- 
sured in Jy, cm, Mpc, and days, respectively. 
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The high radio flux frequently measured in IDV 
sources implies an extremely high brightness temper- 
ature, often many orders of mag nitude above 10 12 K. 
iKellermann fc Paulinv-Tothl (|1969l ) have shown that, as- 



suming the electron distribution follows a single power 
law, the luminosity of the inverse Compton scattered pho- 
tons exceeds that of the synchrotron photons when the 
brightness temperature of the source reaches - 10 12 K. 
Above this threshold, rapid cooling of the relativistic elec- 
trons due to inverse Compton scattering — the "Compton 
catastrophe" — forbids a further increase in the bright- 
ness temperature (see e.g. lKellermannll2002l for a recent 
review of the brightness temperature problem) . The limit- 
ing value is even lower, Tb < 10 11 K, if the magnetic field 
and particle e nergy density of t he source is driven towards 
equipartition (jReadheadl Il994l ) . The observed variability 
in some sources can be interpreted as the result of extrin- 
sic effects, which, at first sight, relaxes the size constraint. 
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For example, the flux variations of PKS 1519—273 and 
PKS 0405—385 are convincingly identified as interstellar 
scintillation. Nevertheless, all realistic models of the scin- 
tillation mechanism impose a new constraint on the size 
and require a brightness temperature of Tr > 10 13 K in 
some cases (jMacquart et al.ll2000t iRickett et al.ll2002h . far 
exceeding the limit imposed by the Compton catastrophe. 

A prevalent feature associated with IDV sources is 
a flat or inverted spectrum (a < 0, with flux F v oc 
v~ a ) at radio-millimeter waveleng ths (e.g., I Gear et al 



1994lKedziora-Chudczer et al.ll200ll ). Optically thick syn 
chrotron emission from power-law electrons rises as j/ 5 / 2 , 
too fast to account for the observed spectra. Optically 
thin synchrotron emission in the scope of the conven- 
tional interpretation of the synchrotron theory has a flux 
F v oc i/ - '' 5-1 )/ 2 , where s is the power-law index of the 
electrons (dNe/dj oc 7~ s ). If a = (s — l)/2 < 0, the num- 
ber density of electrons diverges towards high 7. Imposing 
a high-energy cut-off in the electron spectrum avoids the 
divergence and may account for the commonly observed 
spectr al steepening at optical frequencies, but iMarscher 
( 19771 ) showed that electron spectra with s < 1 would re- 
sult in a high flux between infrared and optical frequencies 
that is not supported by observations. The most common 
interpretation of the flat or inverted spectra is, therefore, 
a superposition of many s ynchrotron spec t ra within an in- 
homogeneous source (e.g. de Bruvn 1976t Marscher 1977 ; 
Blandford fc Konig]|ll979h 

In IKirk fc Tsand (j 20061 ). we discussed a synchrotron 
self-Compton model in which the electron distribution is 
monoenergetic. The lack of low-energy electrons enables 
more GHz photons to emerge from the source, allowing 
a higher brightness temperature to be observed without 
initiating catastrophic cooling. We found that a tempera- 
ture of up to Tb ~ 10 14 K at GHz frequencies is possible 
with only a moderat e Doppler boosting factor of ~ 10. In 
Tsang fc Kirkl (|2007l ), we discussed the parameters of the 
monoenergetic model and showed that the assumption of 
equipartition of energy in the source does not prevent the 
Compton catastrophe. We also showed that an injection 
of highly relativistic electrons or strong acceleration in the 
source cannot produce temperatures much higher than our 
limit due to copious electron-positron pair production. 

In this paper, we examine the spectral properties of 
synchrotron emission from monoenergetic electrons and 
from an electron distribution that is a double power law 
in energy, by comparing the model spectra with the ob- 
servations of S5 0716+714, a BL Lac object that is one of 
the brightes t known IDV sources , as well as a gamma- 



ray blazar ( Hartman et al. 1999f ). In doing so, we as- 
sume that the dominant targets for inverse Compton 
scattering are produced within the source (SSC model). 
The emission from gamma-ray blazars can also be in- 
terpreted in the context of models in which the tar- 
get photons are created externally (EC model), for ex- 
ample in the broa d line region, the accretio n disk, or 
a molecular torus ( Sokolov fc Marscherll2005f) . However, 
in many sources there is no observational evidence of a 



significant external photon source. This is the case for 
S5 0716+714, where, despite much effort over the past 
three decades, no em ission lines have been detected (e.g., 
Bvchkova et al.ll2006[ ) . Furthermore, XMM- Newton obser- 
vation s of S5 0716+714 in 2004 a nalysed by iFerrero et al 



( 20061 ) and Foschini et a l. (2006) show two spectral com 
poncnts in the 0.5 — lOkeV band, whose variability prop- 
erties appear to favour the SSC interpretation. The re- 
cent extensive simultaneous obser vations of this objec t 
from radio to optical frequencies bv lOstorero et al. (2006), 
together with INTEGRAL pointings at GeV 7-ray ener- 
gies during the same period, provide the best test for our 
model. 

In the following, we present the computation of the sta- 
tionary electron distribution and the resulting synchrotron 
and inverse Compton spectra. The model spectra com- 
puted using th e monoenerget i c elec tron approximation, 
as described in iTsang fc Kirk ( 2007 ). are presented first. 
Although adequate for the radio emission, the monoen- 
ergetic model cannot reproduce the entire spectrum of 
S5 0716+714. We therefore investigate an electron dis- 
tribution that is a double power law in energy — a hard 
low-energy part that softens to a high-energy tail above 
a characteristic energy. In this way, the inverted optically 
thin radio emission is retained and complemented by non- 
thermal synchrotron emission from the high energy tail. In 
section [2 we briefly describe these injection models. The 
resulting stationary electron distribution is calculated in 
section [3] and used for the computation of the synchrotron 
and inverse Compton spectra. In section [5J we compare 
the predictions of these models with the observed spectral 
energy distribution (SED) of the source to S5 0716+714. 
Our findings and some limitations of our approach are 
discussed in section [5] and our conclusions presented in 
section [6] 

2. The model 

The h omogeneous monoenergetic model discusse d previ- 
ously (jKirk fc Tsangl 120061 : ITsang fc Kirkl l2007t ) can be 
completely characterised by the Doppler boosting factor 
V = l/[r(l — flcos-d)] (cf3 is the source speed with re- 
spect to the rest frame of the host galaxy, ■& the an- 
gle between the velocity and the line of sight, and L = 
(1 — /3 2 ) -1 / 2 ), the redshift of the host galaxy z, and four 
source parameters, the electron number density N c , the 
magnetic field strength B, the linear size of the source 
R, and the electron Lorentz factor measured in the rest 
frame of the source 7. For the purpose of comparison 
with observations, these can be transformed into a dif- 
ferent set of pa rameters. Deta il s of the transformation 
can be found in Kirk fc Tsand (l2006h . in which iV c , B, 
and 7 are replaced by the characteristic frequency of syn- 
chrotron emission, v s = 7 2 ^o, where Vq — 3eB / (Airmc), 
the Comptonisation parameter £, which is the ratio of 
the luminosity of each successive generation of inverse 
Compton scattered photons to the luminosity of the pre- 
vious generation: £ = 47 2 tt/3, (where tt = N c Ra^ is the 
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Thomson optical depth), and the optical depth r s to s 
chrotron self-absorption at the observing frequency. 1 
size of the source, R, can be constrained, for example, 
applying causality arguments to the variation time, At 
the source: R < cAffl/(l + z). 

We present in Section [4] the model spectra from i 
noenergetic electrons that show good agreement with 
observations of S5 0716+714 at radio frequencies. The 
tical data can be fitted by this model if a high-ene 
power-law "tail" is added. To do this, we consider 
injection spectrum of the form Q(j) oc ('j/^p)~ s 
Ymin < 7 < 7max, where the power-law index s eqr 
si for 7 < 7 P , and S2 for 7 > 7 P (Fig. [I]). The elec- 
tron number density at a given time is proportional to 

7max for S < !) K 7m7n for S > ^ and X m 7max for 

s = 1. In the high-energy branch of the injection spec- 
trum, for 7 > 7 P , we require that S2 > 1, so that electron 
number density congregates towards 7 P . In the low-energy 
branch, 7 < 7 P , the electrons congregate at 7 P if si < 1. 
But we also require that the opacity to synchrotron self- 
absorption is dominated by electrons with 7 — 7 P , which 
is achieved by demanding si < 1/3. Under the conditions 
s\ < 1/3 and S2 > 1, the low- frequency synchrotron spec- 
trum is well-approximated by that of monoenergetic elec- 
trons with Lorentz factor 7 p . 

The electron injection spectrum cuts off at 7 m ; n to- 
wards low energy and at 7 max towards high energy. The 
exact value of 7 m i n is unimportant, since, as explained 
above, synchrotron emission and opacity are dominated by 
electrons with 7 = 7 P in the low-energy part of the injec- 
tion spectrum, where 7 max determines the high frequency 
cut-off in the synchrotron spectrum, at v max = 7max t/ 0j 
and the highest photon energy achievable through inverse 
Compton scattering in the Klein-Nishina limit, which 
equals jm^mc 2 . 

To summarise, the injection spectrum has the form 



Q(7) = Qo 



(l/jp) 31 , 7min < 7 < 7p 
max 



(2) 



where Qo is the electron injection rate per unit volume 
per unit 7 at 7 = 7 p . 

3. Stationary solution 

The shape of the synchrotron spectrum is determined by 
the stationary electron-energy distribution. Electrons in- 
jected into the source according to Eq. © are subject to 
radiative cooling while in the source and evacuate this 
zone on a time-scale close to the light crossing time, 
t csc ~ R/c. The evolution of th e electron spectru m is gov- 
erned by the kinetic equation ( Kardashevlll962 ): 



dn c ( 7 \ d . n c 

Ot \7p/ °1 tcsc 



(3) 



where, for simplicity, we denote the differential electron 
number density (dNc/d'y) by n e . The second term on the 
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Fig. 1. Schematic representation of the electron injection 
spectrum and the stationary differential number density 
as a function of 7. The heights of the spectra have been 
adjusted for easy comparison and are not to scale. The 
solid line shows the double power-law injection spectrum 
with power-law index s\ for 7 < 7 P , and S2 for 7 > 7 P - The 

dashed line shows the case where 7 coo i = 7^1 > 7 P - The 
differential electron number density n c oc 7~ 51 in regions 
A and B, n c oc 7 _S2 in region C, and n c oc ^~( S2+1 ) in 
region D. The dotted line shows the case where 7 CO oi = 
7cooi < 70; with n c oc 7~ Sl in region A, n c oc 7~ 2 in 
region B, and n c oc j^( s ^+ 1 ) in regions C and D. 

right hand side of Eq. ([3]) is the rate of change of the 
electron Lorentz factor due to radiative losses. This term 
is the sum of the rates for synchrotron losses and for losses 
from inverse Compton scattering: 



Ttotal = 7s + Tic 
where 



7s 



4:(t t Ub 2 
3mc 



(4) 



(5) 



The third term is the rate at which electrons escape from 
the source. 

In the stationary state, Eq. Q can be solved analyti- 
cally: 



1 



Q(Y) 

T total 



/i(7) 

with the integrating factor 



/i(7")d 7 " 



/j( 7 ) = 7 exp 



(7'totalW) d 7 ' 



(6) 



(7) 



However, Eq. ([3]) is only a rough description of a 
source, for example, because of the crude treatment of 
particle escape involved in setting t csc — R/c. Therefore, 
rather than use Eq. ^ , we choose to use an approximate 
solution that more clearly demonstrates the effects that 
cooling and the evacuation of electrons from the emission 
region have on the electron energy distribution. 

We first identify the Lorentz factor, 7 CO oi, which de- 
termines the electron energy at which radiative cooling 



4 



O. Tsang and J. G. Kirk: SSC model for S5 0716+714 



region: 

7total 



dominates losses due to particles escaping the emission n c oc 7~ s , since i csc is independent of particle energy. 

On the other hand, when synchrotron losses dominates, 
^ such that t csc > tcooi, the stationary solution is n c oc 7~ 2 

= . (8) for 7 < 7 p , and n e oc 7~( s+1 ) for 7 > 7 p . For the com- 

7=7cooi ^ osc putation of the low frequency synchrotron emission, the 

distribution can be approximated by a monoenergetic one 
at 7cooi in the first case and 7 P in the second. 

The iterative procedure used to find 7 CO oi is as follows: 
The loss rate is defined as 



7 



In principle, 7 coo i can be evaluated only if the entire elec- 
tron distribution is already known, since 71c depends on 
the spectrum and intensity of emitted radiation. However, 
in practise, a simple iterative scheme enables it to be found 
rapidly in all the cases we have computed. Assuming it is 
known, solutions of Eq. |3]) that are valid in the limits 
7 "C 7cooi and 7 ^> 7 CO oi are easily found. In the first case, 
cooling is unimportant, and it immediately follows that 



= t csc Q Q (j/jp)' 



for 7 < 7 cool . 



(9) 



In the second, escape is unimportant, and the appropriate 
solution is found by integrating the kinetic equation once: 



7 



Qo [i 'hp) S dV for 7 > 7coo i. (10) 



These solutions intersect close to the point 7 = 7 coo l- Our 
approximation consists in adopting the solution without 
cooling given in Eq. ^ for all Lorentz factors below the 
intersection point and the solution without escape given 
in Eq. (fl"Tj|) for all Lorentz factors above the intersection 
point. 

In addition, we assume and verify a posteriori (see 
Section[5]) that 71c can be approximated by the expression 
for inverse Compton scattering of the synchrotron photons 
in the Thomson regime: 



Tic 



4g T l/ s 2 

3mc ^ 



(11) 



where U s is the energy density of synchrotron photons in 



the source. In this case, 7totai 
approximate solution is 

t esc Q(l) 

-2 r°° A^inr^i 



7 / (7cooii eS c), and our 



7 < O7cooi 



£csc7cooi7 2 L &l'Q{l') a7cooi < 7 



(12) 



where a (~ 1) is determined by requiring the solution (but 
not its first derivative) to be continuous. 

The Lorentz factors 7 CO oi and 7 P give rise to breaks 
in n e , which correspond to the breaks in the synchrotron 
spectrum at v v — 7^0 and v coo \ — 7c OO i^0- Notice that, if 
s < 1 (as in the injection spectrum below 7 P ), n c is approx- 
imately proportional to 7 _1 oc 7 -2 , whereas if s > 1 (as 
in the injection spectrum above 7 P ), n c is approximately 
0C7 -( S +i). 

Two types of stationary spectra result from Eq. (fl~2|) , 
depending on whether the peak of the injection spec- 
trum, 7 P , is below or above 7 coo i- Figure [T] shows the in- 
jection spectrum as a solid line, the stationary spectra 
where 7 P > 7 coo i as a dotted line and where 7 P < 7 CO oi 
as a dashed line. When electrons are predominantly re- 
moved from a certain energy range by leaving the source 
(t csc < icooi), the spectrum retains its original shape, 



Ttotal = Z 7- 

Amc 



(13) 



Then, starting with (5 = 0, 7 coo i is evaluated from Eq. {8} 
and, using the electron distribution given by (|12p . U s is 
evaluated as described in Sect. 13.11 The value of 5 is read- 
justed to 6 — U s /Ub and the cycle repeated until succes- 
sive values differ by less than 1%. In the examples dis- 
cussed in this paper, convergence was achieved after two 
iterations. Because the change in 7 coo i between iterations 
was only roughly a factor of 2, the final emission spectrum 
was close to that found using 6 = 0. 

3.1. Synchrotron and inverse Compton emission 

The synchrotron specific intensity, following straightfor- 
wardly from the radiative transport equation, is 

1^ =S V {1- exp(-r s )] (14) 

where the optical depth to synchrotron radiatio n is r s = 
g v ■ R , and a v is the absorption coefficient (e.g.. lLongair 
19921 Chapter 18) 



3%/3 



or mc z/l sin < 



16 Of hv 



v 



C Vn <(^h ,15) 

where af is the fine structure constant, <fi t ne angle be- 
tween the magnetic field and the direction of the emit- 
ted radiation, x = v/{^ 2 vq), F{x) =x J K 5 / 3 (t)dt, and 
K 5 / 3 is the modified Bessel function of order 5/3. The 
source function is 



S u = —2mv 2 



/ 7 7 —^(a:K(7)d7 



(16) 



In the monoenergetic approximation, the source function 
simplifies to 



S v = mvl 



-f 5 F(x) 
K 5 / 3 (x). 



(17) 



Equation (fT4j) is integrated over frequency and angle to 
give the energy density of synchrotron photons in the 
source 



U s 



(18) 



where C is a geo metrical factor that is shown in 
Tsang fc Kirkl (|2007l) to be C = 2/3. 
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The synchrotron photons are repeatedly scattered by 
the energetic electrons to higher energies. Denoting by 
i the number of times a photon is scattered, the rate 
of scattering the (i — l) th generation of photons into 
the frequency interval dvj by a single electron (see e.g., 
Georganopoulos et al.ll200ll Eq. (4)) is 



/ dn ph 

where 
Nu^ = 



3cttc 



sp 



4i><. 



4tt CW 



(19) 



(20) 



c hvi-i 

is the number density of the target photons, and I Vi _ x the 
specific intensity of the (i — l) th generation of photons. 

The first generation of scattered photons is produced di- 

(s) 

rectly from the synchrotro n photons: i = 1, 1 Vq — Ih ■ 
Rvbicki fc Lightmanl (| 19791 Chapter 7) assumed that scat- 
tering in the Thomson regime is isotropic in the rest 
frame of the electron, and obtained f(y) « /iso(y) = 
2(1 — y)/3. Here, we include the Klein-Nishina effects (e.g., 
Georganopoulos et al.ll200ll ). in which case 

(4 ei _ l7 y) 2 (l-y)" 



m = 



2y\ny + y+l-2y' 
xP(l/4 7 2 ,l,y), 



2(l+4 ei _i7y) 



y 



(21) 
(22) 



1 1 - uli) 

are the energy of the target photons 



4e 4 -i 7 2 ( 
where e,_i and ti 
and scattered photon, respectively, in units of me 2 , and 
P(l/4 7 2 ,l,y) = 1 for l/4 7 2 < y < 1, and zero otherwise. 

Assuming a spherical source, the rate of scattering 
photons with energy hvi-i to energy hvi, in the observer's 
frame, from a homogeneous distribution of electrons with 
differential number density n c can be found by integrating 
over the electron energy distribution, 



/ dn ph 



V dtdi>i 



(23) 



sp 



Note that R (the linear size of the source) is divided by 2 
to obtain the source radius. 

The specific intensity of the i th generation photons 
is then the scattering rate of the electron distribution in 
Eq. (|2"3")l integrated over all target photon frequency, 

j(C) ( dE 



dtdvidr 2 d£l 

dn ph 



(24) 



/n v d * d ^ ) 4tt(P/2) 2 ' 

and for a general electron distribution n e , it can be written 

as 



/(f) = \Ro^ 



' d 7 



o r 



1ZA 



For a monoenergetic electron distribution, Eq. 
simplified to 

■J 1 JO v i-\ 



(25) 
5| can be 

(26) 



Equations ([26]) and (f25|) are integrated numerically for 
monoenergetic electrons and for an electron distribution 
given by Eq. |T! 



4. The BL Lac object S5 0716+714 

Observations of S5 0716+714 have shown that the 
source exhibits intr a-day variability in the radio and 
optic al bands (e.g. Ghisellini et al.l Il997 : iRaiteri et al 
I2OO3I ) . Correlation between radio (at 5 GHz) and opti- 



cal (at 650 nm) variability suggest that scintillation, a 
process that is not effective at high radio and opti- 
cal frequencies, d oes not play a large part in the ob - 
served variability (jQuirrenbach et al ll99lllWagneill200lh . 
More recent multi-frequency studies of S5 0716+714 (e.g. 



Ostorero et al.l 120061 : Agudo et al. 20061 ) have obtained 



simultaneous measurements from radio to optical fre- 
quencies during the INTEGRAL pointing period, and 
the non-detection of the source by INTEGRAL has pro- 
vided upper limits at X-ray frequencies. Flux variations 
were detected at 32 and 37 GHz over a period of At = 
4.1 days (in November 2003) when the two- frequency mea- 
surements overlap. Since int erstellar scintillation is inef- 
fective at these frequencies, lOstorero et al. ( 20061 ) con- 
clude that the variability is intrinsic. Assuming Hq = 
70kmsec- 1 Mpc~ 1 , with fl x = 0.7, fl M = 0.3, and Q k = 0, 
and a redshift z > 0.3 based on the non-detection of a host 



galaxy (e.g. lWagner et al.lll996h . they derive a variability 
bri ghtness temperatu re of T V ar > (2.1 ± 0.1) x 10 14 K. 

I Bach et al.l ( 2005 ) analysed the data set of VLBI im- 
ages of 11 jet components of S5 0716+714 at 4.9 GHz, 
8.4 GHz, 15.3 GHz, and 22.2 GHz, observed between 1992 
and 2001. Assuming that all the jet components move with 
the same speed along the jet (i.e. all components have the 
same bulk Lorentz factor), they propose that the observed 
wide range (from 5.5c to 16.1c) of apparent component 
speeds is due to variations in the viewing angle and limit 
the Lorentz factor and the viewing angle of the VLBI jet to 
r > 15 and 9 < 2°, respectively. Under these conditions, 
the range of Dop pler factors wou l d be T > ss 20 — 30. 

According to lOstorero et al. I (120061) . observations of 
S5 0716+714 between 5 GHz and 32 GHz can be fitted 
with spectral indices «5_32 of —0.3 and —0.5 at two dif- 
ferent epochs. They suggest that the radio observations 
can be interpreted as optically thick synchrotron emission 
from an inhomogeneous source, and the spectral break at 
v w 10 13 Hz would correspond to the self-absorption fre- 
quency t'abs- In the near-infrared to optical band, obser- 
vation s from 2001—2004, reported by IHagen-Thorn et al 



(|20061 ) suggest that the spectral energy distribution be- 
tween the frequencies = 1-38 x 10 14 Hz and v-q = 
6.81 x 10 14 Hz can be fitted by the power law F v oc v~ 112 . 

Here, we apply the two homogeneous models described 
in section [2] the monoenergetic one, which successfully 
models the hard radio spectrum with relatively few free 
parameters, and the one with double power-law injection 
(and, consequently more free parameters), which also en- 
ables the high-energy emission to be modelled. We adopt 
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a value of z = 0.3 for the red-shift of S5 0716+7 
and a linear size inferred by the variability time scak 
At = 4.1 days, so that R = cAffi/(l + z). The values 
find for the free parameters of the models and for sev* 
parameters derived from these are given in Table [TJ 1 
spectra predicted by the two models are shown in Fig! [] 
and [3j and are discussed separately in the next two s 
tions. 

Identifying the creation of our homogeneous sou 
with an event that leads to the ejection of an individ 
blob observed with VLBI, we estimate the minimum 
power implied by each set of parameters by multiply 
the total energy content of the source by the average raxe 
at which blobs are ejected from the core. In the co-moving 
frame of the source, for monoenergetic electrons, the total 
energy content is 



•^blob — 



D 1 



+ N' cl mc 2 R 



(27) 



and, for power-law electrons, 



^blob - ( 



n' c (j)jmc 2 d'y I R 



-,/3 



(28) 



Transforming to the rest frame of the galaxy, -Ebiob 
T^biob' so ^at the jet power 



( 1 



blob 



\At c 



(29) 



where Ai C j is the average time difference between the 
ejection of successive blobs. Linear fits of the change in 
the position of the 11 jet components (jBach et al. 2005) 
suggest that a new component is ejected from the core 
every 0.1 — 1.8 years. The lower limits to the jet power 
of each model given in Table [JJ were computed using 
At c ; = 1.8 years. 



4.1. Monoenergetic electrons 

In the monoenergetic model, the spectrum is specified by 
four parameters, 2/ a b s , z/ p , V, and £, as well as z and At 
(which are kept fixed for all models). The self-absorption 
frequency f a b s is determined by the first spectral break at 
~ 4 GHz, and v v corresponds to the spectral cut-off. In 
Fig. [21 we compare two models in which one has a cut-off 
at the spectral break at ~ 10 115 Hz, and the other cuts off 
just before reaching the optical point. The Doppler factor 
T> affects the level of the observed flux both by determining 
the linear size of the source and determining the amount 
of boosting the flux receives. The Comptonisation param- 
eter £ determines the ratio of the synchrotron flux to the 
inverse Compton flux, as well as the value of 7 P . Therefore, 
after specifying and u p , £ must be adjusted to com- 
pensate for its effect on the level of the observed flux and 
to ensure the inverse Compton spectra do not exceed the 
INTEGRAL upper limits, while minimising T>. 

Figure [2] shows the measurements and variable ranges 
obtained from the simultaneous multi-frequency obser- 
vation of S5 0716+714 from the study conducted by 




12.5 15 17.5 
Log 10 v [Hz] 



22.5 



25 



Fig. 2. Spectral energy distribution of S5 0716 +714 



Multi- frequency simultaneous data from lOstorero et al 



(2006) are shown as black symbols. Black dots show data 
points, variation ranges are shown by a vertical bar be- 
tween two symbols, and downward arrows show upper 
limits. Values of the parameters are shown in Table QJ 
The model spectra are computed from a distribution of 
monoenergetic electrons and are shown with solid and 
dashed lines. The dashed line shows the strong-magnetic- 
field model where the parameters are chosen such that it 
goes through the data points at optical frequency, whereas 
the solid line shows weak-magnetic-field model where the 
parameters are chosen to mimic the spectral turning at 
10 115 Hz. The values of the parameters are shown in 
Table [TJ 



Ostorero et al 



( 2006h . Also shown are the spectra pre- 
dicted by the model assuming electrons are monoener- 
getic. The Doppler boosting factor is V = 55 in both 
models. The weak-magnetic-field model has a synchrotron 
self-absorption frequency of f a b s = 3.9 GHz with a syn- 
chrotron spectrum that peaks at v v — 300 GHz (the val- 
ues of other parameters are shown in Table [TJ. This model 
gives a brightness temperature of T-q — 3.9 x 10 12 K at 
v ohs = 32 GHz (T B = c 2 F v /(2k B v 2 e 2 ), where k B is the 
Boltzmann constant and 9^ the angular diameter of the 
source). The synchrotron spectrum shows good agreement 
with the data points at radio frequencies. The first-order 
inverse- Compton spectrum gives emission from optical 
to soft X-ray frequencies and the second-order spectrum 
gives gamma-ray emission up to ~ 40 MeV, while emis- 
sion from higher orders scattering is negligible due to the 
Klein- Nishina effect. 

The strong-magnetic-field model has f a b s = 3.9 GHz, 
with its synchrotron peak at v v = 55 x 10 12 Hz (the values 
of other parameters can be found in Table [TJ and gives a 
brightness temperature of T B — 3.7 x 10 12 K at the ob- 
serving frequency z/ D b s = 32 GHz. The synchrotron spec- 
trum extends up to optical frequencies, and gives a reason- 
able fit at radio frequencies up to ~ 10 115 Hz. The first- 
order inverse- Compton spectrum gives X-ray emission, 
the second-order inverse-Compton spectrum is greatly re- 
duced by the Klein-Nishina effects, and very little gamma- 
ray emission is produced. 
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The spectral break at 10 115 Hz is well-fitted by the 
weak-magnetic-field model. We are unable to obtain a set 
of parameters to allow the first inverse Compton spectrum 
to reproduce the optical data. Simple qualitative analysis 
shows that mimicking the optical data points with the 
first inverse Compton spectrum is inconsistent with ob- 
servation. The level of flux that the first inverse-Compton 
spectrum would require in order to account for the optical 
data is much higher than the synchrotron flux (i.e. £ 3> 1), 
so a large 7 would therefore be required, resulting in the 
spectrum extending to frequencies far beyond the optical 
band. The first inverse-Compton spectrum would there- 
fore exceed the INTEGRAL upper limits, and the very 
high X- and 7-ray flux would result in copious electron- 
positron pair production as the 7-ray photons interact 
with the synchrotron photons. 

Alternatively, one can attempt to include the optical 
data in the synchrotron spectrum, as shown by the strong- 
magnetic-field model. The Lorentz factor of this model is 
higher than 7 CO oi) which implies that the particles lose a 
significant portion of their energy by synchrotron radia- 
tion before they vacate the source, so the electron spec- 
trum will evolve into one that is proportional to 7~ 2 . This 
set of parameters therefore violates the monoenergetic as- 
sumption. Furthermore, the predicted spectrum fails to 
account for the spectral break at v ~ 10 115 Hz, and the 
optical flux is very sensitively to the electron Lorentz fac- 
tor. This model is, therefore, inconsistent. It is apparent 
that, in order to reproduce the observed optical emission, 
a power-law component in the electron spectrum at 7 > 7 p 
must be incorporated, which emits synchrotron radiation 
at a frequency above v v . 

4.2. Injection of electrons in a double power law 

The Comptonisation parameter £ in this model must ac- 
count for the 7 dependence of the electron density. It is 
therefore redefined as 



7 "ed7, 



(30) 



which is consistent with the definition used in the monoen- 
ergetic model. However, for the purpose of fitting the flux 
in the INTEGRAL frequency range, ~ 10 18 - 10 20 Hz, 
we introduce a parameter r p , that determines the ratio 
of flux between v v and 7pfp. Inspection of the solid line 
in Fig. [2] shows that synchrotron photons at the spectral 
break v p w 300 GHz are scattered to 7p^ p ~ 10 17 5 Hz. 
Therefore, r p is a parameter that allows us to specify 
the spectral break v v and then to adjust the flux in the 
INTEGRAL frequency range. The normalisation constant 
in the double power-law injection spectrum, Qq is elimi- 
nated in favour of the parameter r p 



4 , f c 



«- c (7)d7, 



(31) 



monoenergetic model and specifying v a \, s and v p , the ra- 
dio and the X-ray flux can be adjusted with £ = r p and 
T>, and v coo \ can be calculated. 

In the low-Doppler-factor model (Figs. [3]and[3]), we 
attempt to m inimise the D o ppler factor of the source. 



Wagner et al.l ( 1996f ) and lOstorero et al 



Accor ding to 

(|2006l ). the variability displayed by S5 0716+714 is in- 
trinsic, and the variation time At = 4.1 days was mea- 
sured at 32 GHz and 37 GHz. Therefore, we require that 
the model spectrum must agree with the data at these 
two frequencies. We are unaware of correlated simultane- 
ous variability measurement at lower frequencies during 
this campaign and, therefore, allow the model spectrum 
to deviate from the data at frequencies below 32 GHz. At 
the expense of having a lower flux than what is observed 
at frequencies below 32 GHz, we find that the minimum 
Doppler factor required is reduced to T> = 30. 

The power-law indices of the injection spectrum used 
to generate the spectrum of the low-Doppler-factor model 
are si = —2 for the low-energy part, such that electrons 
with 7 < 7 p do not contribute significantly to the syn- 
chrotron opacity, and S2 = 2.60, chosen for constructing 
the spectral shape in the infrared to optical band. The rest 
of the parameters are varied while keeping the Doppler 
boosting factor fixed. To find the limiting case, we have 
chosen the self absorption frequency to be v a hs = 32 GHz 
and find that the minimum Doppler factor that can gen- 
erate a high enough level of flux at 32 GHz and beyond 
is I? = 30. The values of the other parameters can be 
found in Table [T] At the observing frequency of 32 GHz, 
the brightness temperature in the frame of the observer is 
Tb = 1.4x 10 13 K. The frequency at which the synchrotron 
spectrum cuts off does not affect the spectral shape at low 
frequencies. However, z/ max = 2?/(l + z) x ^o7max i s con ~ 
strained by the optical data, which imposes a lower limit 
on f max , and the INTEGRAL upper limits, which impose 
an upper limit on v max . The maximum value is shown 
in this model, where ^ max = 10 18 Hz. This equates to 
7max = 5.45 x 10 5 with V = 30 and z = 0.3. 

The low-Doppler-factor model shows that it is pos- 
sible to interpret the observed variability at 32 GHz and 
37 GHz as coming from o ne of the jet c ompo nents with the 
kinematics described by iBach et al.l ( 20051 ). That would 
require the lower-frequency emission to originates from a 
larger region than what is inferred from the observed vari- 
ability at v > 32 GHz. 

In the high-Doppler- factor model (Figs. [3] and 0| we 
assume the emission at all frequencies — including the 
low frequency (< 32 GHz) radio — originates from a sin- 
gle homogeneous source. This is suggested by the correla- 
tion between the v ariability at 5 GHz and 650 nm observed 
in Fe bruary 1990 ( Quirrenbach et al. 199ll : IWagner et aL 
1996). The parameter r p must be kept small, so that the 



with ^0(7) given by Eq. (|12p . In the monoenergetic ap- 
proximation, r p is equivalent to £• Therefore, using the 



inverse Compton spectra are below the INTEGRAL upper 
limits. This is achieved at the expense of a relatively high 
Doppler factor, at T> — 65. The brightness temperature at 
32 GHz is Tb = 2.5 x 10 12 K. This model also shows the 
minimum value of f max , found to be v max = 1.5 x 10 15 Hz. 
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Table 1. Model parameters corresponding to Figs. fj] and [3] 



Mono (dashed) Mono (solid) Power-law (dashed) Power-law (solid) 



Primary Parameters 










z 


0.3 


0.3 


0.3 


0.3 


At (days) 


4.1 


4.1 


4.1 


4.1 


V 


55 


55 


30 


65 


u p (Hz) 


5.5 x 10 13 


3.9 x 10 11 


2.0 x 10 11 


2.7 x 10 11 


v coo \ (Hz) 


1.74 x 10 12 


1.17 x 10 19 


2.07 x 10 16 


2.14 x 10 19 


t'max (Hz) 






1.0 x 10 18 


1.5 x 10 15 


i 


1Q -2.5 


10 -0.75 






Spr n c\ n t*v Parampfpr*! 

^__7 v." ± c\i y -L cu ex J. ± i v." i. *_ J i. o 
















0.92 


0.86 


i? (pc) 


0.15 


0.15 


0.08 


0.17 


6» d (//as) 


32.5 


32.5 


17.7 


38.4 


7p 


691 


800 


244 


696 


O'cool 


123 


4.39 x 10 6 


7.85 x 10 4 


6.19 x 10 6 


'"/max 






5.45 x 10 s 


5.19 x 10 4 


iV c (cm' 3 ) 


0.02 


0.70 


3.15 


0.32 


B (mG) 


648 


3.43 


34.7 


2.65 


5 






0.97 


0.93 


t/ B /(7par 


1.8 x 10 3 


1.0 x 10" 3 


0.04 


8.4 x 10" 4 


Tb (K) 


3.7 x 10 12 


3.9 x 10 12 


1.4 x 10 13 


2.5 x 10 12 




6 x 10~ 9 


2 x 10~ 7 


5 x 10~ 7 


1 x 10" 7 


Pjet (ergs/s) 


1 x 10 45 


3 x 10 43 


5 x 10 42 


4 x 10 43 



Note: 0d is the angular diameter of the source, [/ pa r the energy density of the particles, and Pj Gt is the jet power in the rest 
frame of the host galaxy, as predicted by each model. The compactness of all four models are negligibly small so not included 
in the discussion. 



The high-Doppler-factor model shows that if the emis- 
sion at all frequencies originate from a single source region, 
it must be boosted by a much higher Doppler factor tha n 



proposed for the jet components by iBach et "al 
Even with a viewing angle $ ~ 0, in which case T> 



(120051). 

a 2T, a 



bulk Lorentz factor of T > 33 is required, suggesting either 
that the source was travelling at a m uch higher speed than 
during the observations analysed bv lBach et al.l (|2005| ) or 
that an interpretation of the superluminal motion in the 
VLBI jet that infers a much higher bulk Lorentz factor 
should be applied. One su ch suggestion is a conical jet 
(jGopal-Krishna et al.ll2006l ). Alternatively, the j et compo 



nents may decelerate as they travel down the jet (Marschcr 



1999 : iGeorganopoulos k. Kazanasl 20031 : Ghisellini et al 



2005). The wide range of superluminal velocities shown 



in lBach et alj (|2005l ) would then be a combined result of 
the variations in speed, as well as of the viewing angle. 

5. Discussion 

The spectra of the four models — two monoenergetic and 
two with double power-law electron injection — shown in 
Figs. [2] and [31 show brightness temperatures at 32 GHz 
that are well above the conventional Compton limit. 
However, due to the lack of low-energy electrons, these 
brightness temperatures in fact lie below the threshold of 
the Compton catastrophe (i.e., £ < 1) and are, therefore, 
sustainable by the source. 

An independent limit on Tb is provided by induced 
Compton scattering, when low-energy electrons couple 
with high-frequency photons. The photon occupation 
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Fig. 3. The spectral energy distribution of S5 0716 +714, as represented in Fig. [21 The model spectra, shown as solid and 
dashed lines, are computed from a quasi-monoenergetic electron distribution in the form of Eq. (|12j) . The dashed line 
represents the low-Doppler-factor model where the Doppler boosting factor is minimised, whereas the solid line shows 
the high-Doppler-factor model where the values of the parameters are chosen to account for all radio and optical data 
points. The dashed gri dline shows the p ositio n of 32 GHz. The values of the parameters are shown in TableQ] Historical 
data, as compiled by I Ostorero et all (2006) and shown as grey symbols, at r adio-to-optic a l freq u encies are from 




<uchr et al. (1981), Waltman et al. (1981), Eckart et all (11982ft. Per 



Gear et all (11994). lHales et al l (Il995^ iDoudas et all (Il996l), 



Saikia et all (Il987l ) . iKuehr fc Schmidt! (Il990h. iMoshir et al l |l990h. lHales et all (Il99ll). iKrichbaum et a! 



evl (Il982h IPerlev et all d 1982b . lLawrence et al 



Rengelink et al 



Rile y et al.l (|l999h . ICohen et all (12 002). Ra iteri et all (120 03): 



"(Il997!). iZharig et all (1199 
UV data from iPian fc Treves! (Il993l) . iGhisellini e~al 



dl997l): X-ray data from iBiermann et alj (|l992h . IComastri et all (Il997!). iKubo et a" 
Tagliaferri et all (|2003h . |Pian et all (|2005f i: and ^-rav data from lMcNaron-Brown et al.l (fl99 
Gollmarl lj2nn1 T 



Giommi et al 



Hartman et al 



1999), 



1999), 



number n p h(^) oc I^/is 3 , which implies that the pho- 
ton occupation number is high at and below the peak 
(where v = i^bs) of the synchrotron spectrum. In the 
presence of low-energy electrons, induced Compton scat- 
tering of photons at frequency ^ a b s to frequencies v < 
Vnhs becom es increa s ingly significant as the synchrotron 
flux grows. ISunvaev! (|l97ll ) showed that this process lim- 
its the brightness temperature at a cert ain frequency to 
r R < mc 2 /(fc R T T ) = 5xl0 9 Kfoi-T T ~ l. lSincell fc Krolikl 



1994! ) demonstrated by numerical simulations that rel- 



ativistic induced Compton scattering limits the bright- 
ness temperature of a self-absorbed synchrotron source 
to T B < 2 x lO^-^+S^ 7 ^ 53 K, where ^ GHz is 



mm is the low 
8 For a con- 



the observing frequency in units of GHz, 7 
energy cut-off in the electron spectrum oc 7 
ventional power-law electron spectrum in which 7 m i n = 1 , 
this gives a limit of T B < 2 x 10 11 K at 1 GHz. 

One might suspect that, at such high brightness tem- 
peratures as are predicted by the models shown here, the 
effect of induced Compton scattering should be significant. 
Qualitative arguments reveal the contrary in our model, 
since the low-energy cut-off in the electron spectrum is 
effectively 7 p , and the occupation number (oc F v /v 3 ) of 
photons at frequencies that would couple with electrons 
at 7 P is negligibly small. Alternatively, when using the re- 



sult of ISincell fc Krolikl (|l994T ) . one finds that the limit 
corresponds to T B < 1.3 x lO 13 ^ 7 " (7 P /10 3 ) 3 / 5 K for 
s = s 2 w 2.6 and 7 min = 7 p . 

The spectral brea k at ~ 10 115 Hz was interpreted by 
lOstorero et al. ( 20061 ) as the result of the change in opacity 
of the source. In our interpretation, the break is a result of 
a corresponding spectral break in the electron spectrum; 
below this frequency, the synchrotron spectrum remains 
optically thin. The self-absorption frequency j/ a b s lies at 
a much lower frequency (~ 4 GHz) in our model, which 
implies a weaker magnetic field and/or a lower electron 
density. This interpretation therefore does not require an 
inhomogeneous magnetic field and electron density. 

The spectral index of the high-energy "tail" of the 
synchrotron spectrum beyond the first spectral break de- 
pends on whether 7 P lies below or above 7 CO oi, since this 
affects the final shape of the electron spectrum, as ex- 
plained in Section [3l If the number of electrons leaving 
the energy 7 p mc 2 is dominated by cooling by radiation, 
the synchrotron spectrum continues from F v oc v 1 ^ be- 
tween j/ a bs and t'cooi, to F y oc v^ 1 ! 2 between v coo \ and 
v v , then Fy oc h>~ s ' 2 / 2 between v v and i/ max , and it is 
cut off exponentially beyond i/ max . In this case, the low 
radio-frequency spectrum resembles that of a monoener- 
getic electron distribution of energy j coo imc 2 . If, on the 
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Fig. 4. The spectral energy distribution of S5 0716 +714 
and the model spectra, as represented in Fig. [2j in 
the radio-to-optical band. Top panel shows the low- 
Doppler-factor model where the Doppler boosting factor 
is minimised. Bottom panel shows the high-Doppler-factor 
model where the values of the parameters are chosen to 
account for all radio and optical data points. 



other hand, losses are dominated by electrons evacuating 
the emission zone on a time scale of t esc = R/c 7 the syn- 
chrotron spectrum then continues from F v oc j/ 1 / 3 between 
^ a bs and z/ p , to F v oc ^^( s 2- 1 )/ 2 between v p and f CO oi, then 
F v oc t/~ S2 / 2 between v coo \ and t / maX 7 and it is again cuts 
off exponentially beyond i/ max . The low radio- frequency 
synchrotron spectrum can be approximated in the same 
way as the one from a monoenergetic distribution of elec- 
tron of energy 7 p mc 2 . 

Observations of S5 0716+714 from infrared to opti- 
cal frequencies suggest that the spectral energy distribu- 
tion between the frequencies fx = 1-38 x 10 14 Hz and 

= 6.81 x 10 14 Hz can be well-fitted by the power law 
F v oc j/~ 112 (jHaeen-Thorn et al.ll2006h . Clearly, the top 
panel in Fig. [4] is much too hard at these frequencies. 
However, the spectrum can be softened by lowering the 
cut-off frequency of the synchrotron spectrum ^ max (i.e., 
bottom panel in Fig. [4}. By decreasing v max to approxi- 
mately i/k, the spectrum begins an exponential drop at or 
just before reaching the relevant frequency range and, as 
a result, softens the spectrum, without altering the level 
of flux or the spectral shape at frequencies < u max . 

Figure [2] and the bottom panel in Fig. [4] demonstrate 
that, if the variability of S5 0716+714 is intrinsic and ex- 
tends down to < 32 GHz, the Doppler-boosting factor of 
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Fig. 5. Cooling rates normalised to the synchrotron loss 
rates (dashed lines at 7 = 7 S ). The rate for inverse 
Compton scattering is shown as a dotted line, the total 
loss rate as a thin solid line (blue online) . The approximate 
rate used to compute the electron distribution is shown as 
a thick solid line (black online). The top (bottom) panel 
shows the results of the low(high)-Doppler-factor model 
shown as dashed (solid) lines in Fig. [3J The vertical dot- 
ted lines correspond to 7 = a7 coo i (see Eq. (fT2")l ) on the 
left, and to 7 = 7 max on the right. 



the emission r egion must b e muc h higher than the 20 — 30 
suggested by iBach et al. ( 2005 ). One possibility for en- 
abling the Doppler factor to fall within this range is to 
increase the magnetic field strength and the electron den- 
sity, as shown in Table [TJ This causes the synchrotron self- 
absorption frequency ^ a b s to increase, as shown by the top 
panel in Fig.[4j therefore requiring the assumption that the 
emission at frequency below j/ a b s originates from a larger 
region than inferred from the variability at 32 GHz and 
37 GHz. Requiring the value of v a b s < 32 GHz, we find 
that, to remain below the INTEGRAL upper limits, we 
require a minimum Doppler factor of T> — 30 and a min- 
imum self-absorption frequency of j/ a b s = 32 GHz. Below 
these minima, the model spectrum would either have a 
level of flux below the measured flux at radio frequencies 
or the subsequent first inverse Compton spectrum would 
lie above the INTEGRAL upper limits. 

The models shown here depart from the equipartition 
of magnetic and particle energy, and are dominated by the 
energy of the relativistic electrons (except for the rejected 
model) . Estimating the energy required by the source from 
the host galaxy, Pj ct , we find that the power required from 
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the source by our models is consistent with what is ex- 
pected from a low-en ergy peaked BL Lac object (see e.g., 
Nieppola et"al~ll2006l ). 

An approximation inherent in our method is that, as 
far as their effect on the electron distribution is concerned, 
inverse Compton losses are dominated by single scatter- 
ings off the synchrotron photons and may be treated in 
the Thomson approximation. For each of the models pre- 
sented in Fig.[U we show in Fig. [5] the total rate of radia- 
tive cooling of an electron (thin solid lines, blue online) 
computed using the full Klein- Nishina loss rate for scatter- 
ing on the full output spectrum. This should be compared 
with the thick solid lines (black online) that give the loss 
rates used in computing the electron density according to 
Eq. ([3]). These lines are horizontal and give the quantity 
7totai/7s = 1 + 6 for the converged model solutions. The 
range over which cooling is important in our approximate 
solution lies between the vertical lines at 7 = <27 coo i and 
7max in this figure. (In the lower panel, these lines lie close 
together.) The maximum deviation in this range is roughly 
20% and occurs at 7 = 7 max for the model with a relatively 
low Doppler factor (plotted as a dashed line in Fig.|4j). The 
deviation for the higher Doppler factor model is less than 
3%. We conclude that this approximation has a negligible 
effect on the electron distribution. The computation of the 
output spectrum from the electron distribution is not af- 
fected by this approximation, since the full Klein-Nishina 
expression for the emissivity is used. 

A testable prediction of our models is correlated vari- 
ability. The hard a w —0.3 radio spectral component be- 
tween 30 GHz and several hundred GHz is interpreted as 
the optically thin emission of electrons whose character- 
istic synchrotron frequency lies at or above hundreds of 
GHz. Consequently, fluctuations in the number of such 
particles should be simultaneously reflected in broad-band 
fluctuations of the specific intensity of radiation at these 
frequencies. We are not aware of studies that investigate 
such an effect. 

The variability of the synchrotron spectrum should 
also be correlated to the variability of the inverse Compton 
spectrum, since the same electrons are respo nsible for both 
emissi ons. Analysis of XMM observation bv lFerrero et al 



(2006) show correlated variability between two X-ray spec- 
tral components of S5 0716+714, which were interpreted 
as the high-energy part of the synchrotron spectrum and 
the low-energy part of the inverse Compton spectrum. 
In the context of our model, the synchrotron component 
is the emission from the high-energy tail of the electron 
spectrum, whereas the IC component is the synchrotron 
emission of the hard low-energy electron spectrum scat- 
tered by the same low energy electrons. The study by 
Ferrero et al. ( 20061 ) therefore does not provide a direct 
test of our model, although their data support an SSC 
interpretation. 

Similarly, the polarisation properties of the 30 - 
100 GHz emission are predicted to be those of the single- 
particle synchrotron emission. One aspect of this is the 
possibility that the intrinsic circular polarisation of the 



source might reach a few percent ( Kirk fc Tsang 20061 ). 
Another concerns the degree of linear polarisation. Power- 
law electron distributions radiating in a source with a 
completely homogeneous magnetic field can theoretically 
reach quite high degrees of linear polarisation (up to 
« 70%, depending weakly on the power-law index). These 
values are not reached in extra-galactic radio sources, 
probably because the magnetic field direction within the 
source is tangled. In general, higher degrees of polarisation 
are found at higher frequency, consistent with the con- 
ventional picture in which the source size decreases with 
increasing frequency. 

In our model, however, the degree of linear polarisa- 
tion in the radio is determined by monoenergetic electrons. 
According to standard synchrotron theory, this tends to 
50% at low frequency (x — vjv s <C 1), rises to 76% at 
x = 1, and tends to 100% at x ^> 1. The effect of a tan- 
gled field within the source reduces these values, but is 
independent of observing frequency. Thus, as in conven- 
tional inhomogeneous models, the degree of linear polar- 
isation is predicted to increase with frequency. However, 
more quantitative predictions would require consideration 
of effects, such as internal Faraday rotation, and lie out- 
side the scope of the present paper. 

Finally, we note that the hard electron injection spec- 
trum that we have adopted (n c tx 7 2 ) corresponds to 
the low-energy (7771c 2 < k^T) part of a relativistic ther- 
mal distribution. According to Table 1, the corresponding 
electron temperature would lie at around 10 12 K. In an 
electron-ion plasma, this corresponds to the temperature 
of shocked gas behind a mildly relativistic shock front, if 
one assumes that the electrons and ions equilibrate to a 
common temperature. Recent P.I.C. simulati ons suggest 
that this assumption may indeed be justified (jSpitkovskv 
2007h . 



6. Conclusion 

Using the specific case of S5 0716+714 as an example, 
we confirm that it is possible to produce high brightness 
temperatures at GHz frequencies in compact radio sources 
without the onset of catastrophic cooling, provided that 
the radiating particles have a distribution that is suffi- 
ciently hard below a characteristic. In addition, we show 
qualitatively that induced Compton scattering is insignif- 
icant in sources with a low-energy electron cut-off despite 
the high brightness temperature, the underlying reason 
being the low occupation number of the photons that can 
couple with the electrons at the cut-off energy. 

The model where an electron distribution that is a 
double power law in energy, peaking at 7 P , is injected into 
the source offers more flexibility at higher frequencies in 
the synchrotron spectrum (from infrared to optical) at the 
expense of more free parameters, compared to either mo- 
noenergetic or single power-law distributions. These pa- 
rameters should be constrained by simultaneous observa- 
tions due to the highly variable nature of IDV sources. In 
the case of S5 0716+714 where such data is available, the 
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spectral break at about 230 GHz determines the value of 
7 p , the optical data at 5 x 10 14 Hz gives the lower limit of 
^maxj a $ well as constraining the spectral index S2, and the 
INTEGRAL upper limits give the upper limit of ^ max and 
also constrain the value of r p , which in turn determines 
the electron density. 

The example of S5 0716+714 illustrates several impor- 
tant spectral properties of an electron distribution with a 
low-energy cut-off, as described in the previous sections. 
The most noticeable feature is the hard, inverted, op- 
tically thin synchrotron spectrum, spanning a wide fre- 
quency range, which is a prevalent fe ature in compact 
radio sources at radio frequenc ies (e.g-. lGear et al.lll994t 
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200 lh 

spectral breaks at i/ p = j^Uq, i^cooi 
exponential cut-off at ^ max = 7^ lax i / c 
fore, allows a simple homogeneous source to reproduce the 
common features shown by many IDV sources. 



Other features are the 

= Tcool^O: and thc 

This model, there- 
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